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We report on our recent measurements of the Cabibbo-Kobayashi-Maskawa CP-violating phase 7 and of related CP- 
asymmetries and branching fraction ratios. The measurements have been performed on samples of up to 465 million 
BB pairs collected by the BABAR detector at the SLAC PEP-II asymmetric-energy B factory in the years 1999-2007. 

1. INTRODUCTION 

o : 

Within the Standard Model CP violation arises from a single irreducible complex phase in the Cabibbo-Kobayashi- 
q Maskawa (CKM) quark-flavor- mixing matrix Uand manifests itself as a non-zero area of the Unitarity Triangle (UT). 

OThis triangle depicts, in the complex plane, the relation 1+ Yr d Y/i b + Yf d Yft =0, that follows from the unitarity of V. 

\& , The BABAR experiment has measured precisely, in B — » (cc)X s decays, the angle (5 = arg — Yj d Y/» of the UT, finding 

VtdV tb 

a value that differs significantly from and it. This clearly indicates that the area of the UT is non- vanishing, thus 
^ ' establishing CP-violation. In order to confirm that the CKM mechanism is the correct explanation for CP-violation, 
^ • we need to overconstrain the UT by measuring precisely the other angles (a and 7) and the sides. 

or 

fi: 2. MEASURING 7 WITH B MESON DECAYS IN BABAR 

The angle 7 = arg — Y/ d Y^* can be measured in a theoretically clean way in CP- violating B meson decays to open- 
charm final states, I)WX S . In these decays the interference between the tree amplitudes b — > cus and b —> ucs leads 



X 



A(b->u) 
A(b^c) 



and the relative strong 



O ■ 

iy-j , to observables that depend on the relative weak phase 7, the magnitude ratio tb — 

phase 5b between the two amplitudes, tb and 5b depend on the B decay under investigation; they are difficult to 
CO ' estimate precisely from theory (QCD), but can be extracted directly from data together with 7. 

In order to search for the interference between the b — > cus and b — ► ucs tree amplitudes, we have reconstructed the 
following decays, which lead to final states that can be produced both through b — > c and b — ► u mediated processes: (i) 
charged B -► £>W° /D { * } ° K and B -> D°/TfK*(K* -> K%tt) decays and neutral B° -► D°/TfK*°(K* -> K+w') 
decays, followed by meson decays to final states accessible also to (ii) neutral B° /B° decays to D ± K°tt t . 



In the latter case, interference between the b— >c and the b— >u mediated processes is achieved through B°—B° mixing. 

The measurements presented here exploit, partially or in total, the large data sample (465 x 10 6 BB pairs, assumed 
to be equally divided into B + B~ and B°B ) accumulated by BABAR in the years 1999-2007. Yet, they are still 
statistically limited, as the effects that are being searched for are tiny: 

• The branching fractions of the B meson decays considered here are on the order of 5 x 10~ 4 (B -> D^°K, 
D°K* and DK°n) or 5 x 10~ 5 (B° — > ifx* ). Branching fractions for I)W decays, including secondary 
decays, range between O(10~ 2 ) and O(10~ 4 ); D~ candidates are reconstructed in K~tt + it + (B w 10%). 

• The interference between the b — > c and b — > u mediated B decay amplitudes is predicted to be low, as the ratio 
r B is expected to be around 40% for neutral B° — > D°K*° and DK°ir (due to CKM factors) and around 10% 
for charged B — > £>(*)° }{(.*) decays (due to CKM factors and the additional color-suppression of A(b — > u)). 

The final states of these B meson decays are completely reconstructed, with efficiencies that range between 40% (for 
low- multiplicity, low-background decay modes) and 5% (for high- multiplicity decays). The selection is optimized in 
order to maximise the statistical sensitivity S/ y/S + B, where the number of expected signal (S) and background 
(B) events is estimated from simulated samples and data control samples. 
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3. MEASUREMENTS OF 7 WITH CHARGED B MESON DECAYS 
3.1. Dalitz-plot method 

On a sample of 383 xlO 6 BB pairs we have selected B^D°K, D*°K (D"°^D°j and D°7r°), and D°K*(K*—>KgiT) 
decays, followed by D° decays to the 3-body self-conjugate final states Kgh + h~ (h = n, K). [l[ We have reconstructed 
163 B candidates with D° -> K%KK and 979 B candidates with D° -> K%-k-k. The channel B -> D°K*,D° -> 
KgK + K~ has not been reconstructed due to lack of statistics. Following the approach proposed in [2j, from a fit to 
the Dalitz-plot distribution of the D° daughters we have determined 2D confidence regions for the variables (x±, y±), 
which are related to 7 by x± = tb cos(8b ± 7) and y± = tb sin(<$B ± 7). In the fit we have used a model for the D° 
and D decay amplitudes to Kgh + h~ described as the coherent sum of a non-resonant part and several intermediate 
two-body decays that proceed through known K gh or h + h~ resonances. The model has been determined from large 
(« 5.6 x 10 5 ) and very pure (~ 98%) control samples of D mesons produced in D*^Dir decays. 

The results for x and y are summarized in Table [T] They are consistent with and have similar precision to 
those obtained by the Belle experiment. (s| From the (x±,y±) confidence regions we determine, using a frequentist 



Table I: Measurements of x± and y± obtained with the Dalitz and GLW analyses of B — > ' K 



parameter 


B -> D°K 


B -> D*°K 


B -> D°K* 


x+ (Dalitz) 
x+ (GLW) 
y+ (Dalitz) 
x- (Dalitz) 
x- (GLW) 
y- (Dalitz) 


-0.067 ± 0.043 ± 0.014 ± 0.011 

-0.09 ±0.05 ±0.02 

-0.015 ± 0.055 ± 0.006 ± 0.008 

0.090 ± 0.043 ± 0.015 ± 0.011 

0.10 ±0.05 ±0.03 

0.053 ± 0.056 ± 0.007 ± 0.015 


0.137 ± 0.068 ± 0.014 ± 0.005 

0.09 ±0.07 ±0.02 

0.080 ± 0.102 ± 0.010 ± 0.012 
-0.111 ± 0.069 ± 0.014 ± 0.004 
-0.02 ±0.06 ±0.02 
-0.051 ± 0.080 ± 0.009 ± 0.010 


-0.113 ± 0.107 ± 0.028 ± 0.018 

0.125 ± 0.139 ± 0.051 ± 0.010 
0.115 ± 0.138 ± 0.039 ± 0.014 

0.226 ± 0.142 ± 0.058 ± 0.011 



procedure, la confidence intervals for 7, tb and 5b (Fig. Q}. We obtain 7 mod 180° = (76^24)°- The total error is 
dominated by the statistical contribution: the experimental and Dalitz-model-related systematic uncertainties amount 
to 5° each. We find values of tb around 0.1, confirming that interference is low in these channels: r^° K = 0.086±0.035; 
r°'" K = 0.135±0.051; kr^" K ' = 0.163±g;?g| (fe=0.9±0.1 takes into account the K* finite width). The small values 
of tb favored by our data are responsible - since cr-y~cr Xi y/rB - for the larger uncertainty on 7 when compared to the 
analogous Belle measurement, which favours values of tb about twice higher than ours. We also measure the strong 
phases (modulo 180°): <5g° K =(lQ9i^)°; 6§*° K =(-63±H)°; 6%° ' K * =(104±g)° . A 3d evidence of direct CP violation 
is found when comparing y+) to {x-, yJ) (they are equal in absence of CPV) in the three B decay channels. 




Figure 1: Two left-most plots: 1-CL as a function of 7 and tb for B —+ D°K, D*°K and D° ' K* as obtained with the Dalitz 
analysis of D° — > K%h + h~ . Two right-most plots: GLW observables measured by BABAR and other experiments. 
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3.2. Gronau-London-Wyler ("GLW") method 

We have reconstructed, on a sample of 383xl0 6 BB pairs, B^D°K and D*°K (D"°^D°j and L>°7r°) decays 
with neutral D mesons decaying to CP-even (K + K~, ir + ir~) and CP-odd (Kgir , Kg(f> and KgU>) eigenstates. 
We have also reconstructed B—*D(*)°K, D°^K~tt + and B-^D^^tt decays as normalization channels. We have 
identified almost 500 B~^D® CP K Q| and 500 B^D*£ P K Q decays. From the observed B ± yields we have de- 
termined the "GLW" observables Q, iZ&L= r ^ B ^ D cll K ~)+ T ( B+ ^ ( £± K+ ) = i +r W 2 ±2r ( * ) cos 7 cos4* ) and 

CP± (r(B--*DM«K-)+r(B+^D ( )0 K+))/2 B B I B 

A^l,=^— D ?^n K ^ D ?5n K ^ = ±2r^ sin 7 sin Jo V-R^p-u- Here we use the notation rl* =rS'' l ° ,i ' and 

6g^=Sg ( } ° K . The results, summarized in Table HT1 and Fig. [I] constitute the most precise measurements of the CP 
asymmetries and branching fraction ratios in these decay channels. A 2.8c hint of direct CPV is seen in B — > D^ p ,K 
decays, consistent with the results from the B~^D°K, D°~^K s h + h~ Dalitz analysis and with recent results of the 

GLW analysis of B^D^ p+ K decays performed by CDF. From the relation r p = CP+ — 2 CF ~ — we have obtained 
loose determinations of r { * ] : (r§° K ) 2 = 0.05±0.07±0.03, (rg '° K f = 0.22±0.09±0.03. We have also determined the 
variables x± by exploiting the relations X j_ = Rcp +( 1 T A cp+)-Rcp-(i : fAcp-) . ^ e resu ^ Sj consistent with and similar in 

accuracy to those from the Dalitz analysis, are summarized in Table |TJ Due to a discrete 8-fold ambiguity on 7 and 

(*) 

the large uncertainty on r B , the GLW results alone do not constrain 7, but are expected to improve the accuracy 
<7 7 of the Dalitz results presented in the previous section by a few degrees when combined with them. 



Table II: Results of the B -> D°K and B -» D*°K GLW analyses 



B decay 




pW 

rx CP+ 


M 

A CP- 


pW 


D°K 
D*°K 


0.27 ± 0.09 ± 0.04 
-0.11 ±0.09 ±0.01 


1.06 ±0.10 ±0.05 
1.31 ±0.13 ±0.04 


-0.09 ±0.09 ±0.02 
0.06 ±0.10 ±0.02 


1.03 ±0.10 ±0.05 
1.10 ±0.12 ±0.04 



4. MEASUREMENTS OF 7 WITH NEUTRAL B MESON DECAYS 

4.1. B° -> D°K*°: Atwood-Dunietz-Soni ("ADS") and Dalitz-plot methods 

In the "ADS" analysis § we reconstruct D° decays to doubly-Cabibbo-suppressed final states: K + n~ , K + n~Tr° 
and K + tt~ 'n + . The Cabibbo-allowed charge- conjugate final states are used as normalization and control sample. 
In principle 7 could be determined from the branching fraction ratios (Rads) and the charge asymmetries {Aads) 
of these decays. 0] In practice we do not have enough statistics to measure Aads and put useful constraints on 7; 
on the other hand by measuring Rads we can infer a significant constraint on rs=r B K °. We have selected, on a 
sample of 465 x 10 6 BB pairs, 24 signal candidates summed over the 3 D° decay channels, and put 95% probability 
bayesian limits on Rads- Rads < °- 244 ; r ads < °- 181 ; r ad¥ < °- 391 - We have then combined those 3 results 
and determined bayesian probability regions for r$ by adding external information (from c- and B-factories) on the 
jD^ decay amplitudes. At 95% probability rs lies between 7 and 41% (Fig. [2]). 

In the Dalitz analysis P we reconstruct D^ 3 decays to K s tt + tt~ . Using 371 x 10 6 BB pairs we have selected 39 
signal candidates. From the fit to the Dalitz-plot distribution of the D° daughters, using the same model for the 
decay amplitudes as in Sec. 13. H we have determined 2D bayesian probability regions for (rs, 7) (Fig. [2]). 
After combination with the rs likelihood from the ADS method we obtain, at 68% probability, 7 mod 180°=(162±56)°. 



-ftO 



4.2. 2/3 + 7 from a time-dependent Dalitz-plot analysis of B°/B 



D ± K°7T* 



On a sample of 347 x 10 BB pairs we have searched for neutral B meson decays to D ± Kgit^ L . [10| We retain 



only those events where the other B meson (Stag) decays to a final state which allows to determine its flavor 
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Figure 2: Left: likelihood function for rs from the ADS analysis of B°— >D°A"*°, including 68% (dark-shaded) and 95% (light- 
shaded) probability regions. Center: 7 vs rs 68% and 95% probability regions obtained from the B°—>D° K*° Dalitz analysis. 
Right: distribution of the fitted values of 2/3 + 7 from selected B° — > DK°tv decays for different hypotheses on the value of r. 



(either B° or B ) and the proper time difference At between the two B decays, inferred from the separation of 
the B decay vertices and momenta. The effective efficiency of the flavor-tagging algorithm is around 30%. A 
fit to the Dalitz-plot distribution of DK®ir as a function of At yields 2/3+7 through the relations: T(x, At, £, 77) = 

' Atl K{l-^C(x)co S (Am d At) HS v (S)sin(A md At)}, S v (x) = ^M^™?™?™-^ , 
tb is the B° lifetime, Am the B° — 5° mixing frequency, x the position in the DKtt Dalitz 



X 



4tb 

A c (S) 2 -A M (x) 2 



C(x) 

plot; £= — 1(1) if the flavor of the B tag is B°(B ), and 77= + 1(— 1) if the final state contains a D + (D~). 

We have selected 558 signal candidates and performed the time-dependent Dalitz-plot fit. For the direct B decay 
amplitude we use a model which incorporates all known intermediate resonances and float most of the magnitudes 
(A) and phases (</>) in the fit. The ratio r between the 6— >w (A u ) and b^c (A c ) amplitudes is fixed. A scan of 2/3+7 
as a function of r is shown in Fig. [2j For the final 2/3 + 7 result we assume r = 0.3 ±0.1. Using the world average for 
j3 from B decays to charmonium we obtain 7 mod 180° = (40 ± 57)°, in agreement with the other determinations. 



5. CONCLUSION 

Many new results on 7 have been obtained by BABAR during 2008. They are in good agreement with each other 
and favor 7«72°, dominated by the B—*D(*)°K*, D°~^Kgh + h~ Dalitz analysis and consistent with expectations 
from CKM fits. Thanks to the combination of several strategies, the uncertainty ct 7 is approaching 20°, which was 
not an original goal of the B factories. er 7 is limited by the available statistics; we expect a 20% improvement from 
the collected data that have not been analyzed yet. The interference effects are confirmed to be small both in charged 
and neutral B decays, thus a much larger statistics is needed to reach a precision of a few degrees on 7. 
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